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We find nonperturbative solutions of a nonlocal scalar field equation, with cubic or ex- 
ponential potential on a cosmological background. The former case corresponds to the 
lowest level effective tachyon action of cubic string field theory. While the well known 
Minkowski solution is wildly oscillating, due to Hubble friction its cosmological counter- 
part describes smooth rolling towards the local minimum of the potential. 



1. Introduction 

The search for signatures of nonperturbative quantum gravity in the history of 
the early universe has motivated several authors to consider cosmological models 



inspired by open string field theory (OSFT), 1,2,3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 - 15 ' 16 the p- 
adic string, 14 . 15 - 17 . 18 ' 19 ' 20 ' 21 or other nonlocal effective act ions. 22 ' 23 ' 24 ' 25 ' 26 ' 27 ' 28 ' 29 



In the first category, the typical starting point is a nonlocal scalar equation with 
polynomial potential, 



where □ is the d'Alembertian operator, m and a are constants, and $ = e <E>. 
This type of pseudo-differential operators can be conveniently manipulated with 
the diffusion equation approach, which was developed and employed, in diverse 
numerical or analytical guises, in Refs. io.n, 12, 14,15,20,30,31, 32, 33, 34,35, 36^ j n brief, the 
method can be outlined as follows: 10,20,33 ' 34 (A) Let the scalar field <j>{t) — > (j>(r,t) 
evolve also along an auxiliary direction r and interpret s as a fixed value of the 
extra variable. (B) Find a solution (j>(0,t) of the local system (s = 0). This is the 
initial condition for a system that evolves in r via the diffusion equation 
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(1) 



(□ + a r )$ = o. 



(2) 
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(C) The solution of the diffusion equation is 

*(r,i) = e rD mt) • (3) 

In particular, the effect of the nonlocal operator e qU is a shift of the auxiliary 
variable r: 

e« n $(r,t) = e- qd ^{r,t) = <f>(r-q,t). (4) 

(D) If the initial condition is chosen to be constant almost everywhere, the final 
configuration $(s, t) obtained by diffusion along r is a smooth function which solves 
(exactly or approximately) the original nonlocal system for some s (here we do not 
consider smooth initial conditions 10 ' 33 ). 

In a recent paper, 36 we found a rolling homogeneous solution to (1): 

* = .(«,,) = 2fc|), <„ 

where z = 2-y/r t and 

1 (a,z 2 ) = 2 ( dTT^e-^ (6) 



Jo 

is the lower incomplete gamma function (see Sec. 8.35 of Rcf. 37 ). The constant a > 
will determine the shape of the scalar field profile. First, we derived an approximate 
scaling formula which expresses powers of the incomplete gamma function as a single 
7 with rescaled arguments: 

[$(a,z)]"w$(na,/«z), (7) 

where 

\T(na + 111 1 /( 2 ™") 
[I> + l)] 1/(2a) 

When n = 1 + 1/a, using the properties of 7 one can show that (5) approximately 
(in the sense of (7)) solves (1) with 

1 \ 2 an 2 2 

1 2 r, to = , g — to , (9) 

Ac y r 

and the universe expands as a power-law: 

□ = _ 5 2 _ . (10) 

The solution obeys the diffusion equation (2) and the action of the nonlocal operator 
in (1) is a rescaling: 

e sD $(a,z) = <f>(a,Kz). (11) 

The level of approximation of the solution (5) is determined by the duplication 
formula (7) and is sufficient to construct easily instantonic tachyon solutions of 
supersymmetric string field theory on flat Minkowski spacetime, (n, a) — (3, 1/2). 
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The approximation scheme related with the duplication formula was extensively 
discussed in Ref. 36 . Its domain of validity was studied in the parameter space 
(n, a). The values taken here lie within this domain. 

In this paper we consider the inverse problem of which cosmology one would 
expect for a given scalar field profile. This perspective is sometimes used in infla- 
tionary or big bang cosmology when one requires a particular profile. For instance, 
in several inflationary models the scalar potential is derived from other known 
variables (scalar profile $(t) and Hubble evolution H{t)) by the Hamilton- Jacobi 
equations; if one is interested in runaway potentials, the inverse problem selects 
simple $(i)'s and H(t)'s from which one can extrapolate the field potential with 
those characteristics. A third application of the inverse problem is in the context 
of ekpyrotic/phantom/braneworld dual models, where a duality between physically 
inequivalent models is established through the scale factor or Hubble parameter. 
The scalar potential is consequently extracted either from the Hamilton-Jacobi 
equations or directly from quantities of the dual model (see Ref. 38 and references 
therein). 

Here, the main difference with respect to the cases above is that we also extract 
the Hubble parameter from the diffusion equation for a given solution. The solution 
is given ab initio because the incomplete gamma function is the only profile we 
know which allows one to easily solve the diffusion equation and (approximately) the 
nonlinear scalar equation simultaneously via the duplication formulae of Ref. 36 . The 
diffusion equation determines univocally a family of Hubble parameters, while the 
algebraic and analytic properties of the incomplete gamma fix the scalar equation, 
and hence the potential. The potentials we obtain are rather simple, thus further 
justifying the method a posteriori. The other independent equation of motion (the 
Friedmann equation) may not be the one of standard general relativity and one will 
have to reconstruct it. However, this is not the focus of this work. 

We discuss the dynamics of two cosmological scalar profiles, one with the cubic 
potential of the bosonic OSFT tachyon and the other for a toy model with expo- 
nential potential. In the former case, the tachyon rolls down the minimum driving 
a superaccelerating expansion. We assume the scalar field to be the only matter 
present in the universe and that it sources the cosmological expansion. In general, 
the expansion can also be driven by gravity itself, like in higher- derivative theories. 
This sector is beyond our scope because we do not solve all of Einstein's equations. 

Although it is not shown whether the solutions thus found have a direct physical 
application, they are worth recording as they illustrate a methodology for treating 
nonlocal nonlinear equations analytically on a curved (in particular, Friedmann 
Robertson-Walker) background. Such is the first motivation for this work. More- 
over, they are explicit realizations of the intuitive idea that the cosmological friction 
heavily affects nonlocal dynamics. This is clear for the solution with cubic potential, 
which can be confronted with its unbounded Minkowski analog (the string tachyon): 
the wild oscillations of the Minkowski solution are absent in the cosmological solu- 
tion because, due to Hubble friction, the field reaches the minimum of its potential 
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in the infinite future. 

The present work has been conceived as a follow-up of our line of investiga- 
tion, and as such it is not meant to be self-contained. We refer to other publica- 
tions for extensive discussions on issues which may be relevant to nonlocal diffus- 
ing theories 10 ' 20 ' 33 ' 34 ' 36 and nonlocal cosmological models. 4,10 Anyway, we wish to 
(re)make a cautionary remark. This class of cosmological models are considered to 
be toy simplifications of a more complete scenario. The gravitational sector is not 
taken into account adequately, since it is basically local; so, one is actually looking 
at models where nonlocality has not been retained systematically. 4 ' 10 This problem 
is still unresolved in the literature for this type of nonlocal operators. From the 
point of view of string theory, the issue regards the closed string background, for 
which there is not as much control and consensus as we have for the open string. 
From the point of view of nonlocal toy models, the authors of Ref. 23 consider a 
scenario of pure gravity. Inspired by that, one might try to implement the nonper- 
turbativc diffusion method also into the gravitational sector. Unfortunately, there 
is a technical reason why we cannot do so, namely, that if all sectors were dressed 
with the same nonlocal operators, then the Hubble parameter (friction term in □) 
would depend on r, thus spoiling Eq. (4). A possible way out is to take different 
diffusion times Ti for each sector of the theory, but we will leave that for the future. a 

2. Cubic Potential (Bosonic Tachyon) 

The case (n, a) — (2, 1) is rather interesting because it corresponds to the tachyon 
of bosonic string field theory. 39 ' 40 ' 41 ' 42 For these values, we can estimate the error 
of approximation by computing 

A max = sup A(z) = sup 

z z 

where LHS and RHS are, respectively, the left- and right-hand side of (7). The 
validity of (12) as a robust assessment of the global error has been discussed and 
proven in Refs. 10 > 20 > 33 > 36 . Any solution of the diffusion and scalar equations is easily 
constrained asymptotically to be such (see these papers for details), so there is no 
issue about late- or early-time deviations from an exact solution. If they are small, 
violations of the equations of motion at some intermediate time t = to can be taken 
into account by adding corrective terms to the solution which vanish asymptotically, 



a There is another class of phenomcnological nonlocal models, introduced in Ref. 25 , which is not 
inspired by string theory but by quantum field theory loop effects, where nonlocality is encoded 
in inverse powers of the Laplace-Beltrami operator. By virtue of a field redefinition similar to the 
one used in f(R) models, one can reduce these models to local ones of scalar-tensor type. 26,27 
Prom that point on, the classical analysis is standard. On the other hand, in our case not only do 
wc not yet have an Ansatz for the gravity action, but the nonlocality we specialize in does not 
allow any field redefinition leading to a simple local system (this point was discussed extensively 
in Refs. 10 ' 34 ). 



LHS - RHS 



LHS + RHS 



(12) 
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but in the present case this is unnecessary because the approximate solution $ 
already captures the correct behaviour at all times. 

This can be seen by looking at the global error associated with $, which is small. 
Plugging in the theoretical value k = 2 1 / 4 sa 1.19, the supremum norm of A is 
A max = A(z 1.23) « 1.8%. Although this error is already acceptable considering 
that we have not employed numerical tools, we can improve it by estimating the 
value of k which minimizes the global quantity 



S = inf 



dzA 2 (z) 



(13) 



where z = O(10) is sufficient. We find k„ 
Fig. 1). 



1.18 and 5 1.4% w A max (see 
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Fig. 1. A with the theoretical (thick line) and numerical (thin line) value of re. Here (n, a) = (2, 1). 



Since the metric is singular at t — 0, the solution $(1, z) = 1 — e 2 is meaningful 
only at either negative or positive times. The equation of motion is 

1 



$ + cr($ - $^) = . 



(14) 



The static potential V($) = cr($ 2 / 2 ~ $ 3 / 3 ) has a local 

minimum at V^(O) = 

and a local maximum at V(l) — a/ 6 > 0. Since — V is the potential of the 
OSFT tachyon (at zero truncation level) on Minkowski, we regard $ as a solution 
of the Riemannian problem, i.e., a problem in Euclidean signature on a curved 
background. The cosmological equations of motion in Euclidean signature are simply 
those where the usual Wick rotation t — > it is performed. This interpretation of the 
equation of motion is necessary only if one identifies the field $ with the OSFT 
tachyon, because the cubic potential is inverted with respect to the usual definition 
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in the theory. If one does not wish to do so, then this problem can be considered to 
be a toy model in Minkowski signature. 

Let t < 0. The Hubble parameter is H = a/a = l/(— 3£), corresponding to 
a scale factor a(t) = (— i) -1 / 3 . The Riemannian universe expands (H > 0) in 
superacceleration (H > 0). At t = — oo (a = 0, big bang) the field sits on the 
local maximum and starts rolling down the minimum. However, it will never reach 
it, since the cosmological friction grows with time and the universe will expand 
indefinitely (a,H — > +oo) b (see Fig. 2). In the absence of a constant term in the 
potential, the equation of state evolves from a stiff matter regime (w = p/p <~ 1, 
where p and p are, respectively, the pressure and energy density of the scalar field) 
to a cosmological constant (w <~ —1). 

In the corresponding local system i.e., the one with a and n fixed while sending 
k to zero, if $ is solution of the Euclidean (in this case, Riemannian) problem, 
then 'I' = 1 — $ is solution to the Lorentzian equation of motion. When nonlocal 
effects are taken into account, this is true only up to a term cx (\I> — vp), which 
vanishes in the local limit. (Nonlocal terms of this form were briefly discussed in 
Rcf. 10 .) Qualitatively, the dynamics of the nonlocal Lorentzian model with this 
"evanescent" term is the same as the one without it. For t < 0, rolls down 
from the local maximum V(0) = towards the local minimum at ^> — 1, driving 
superacceleration of an expanding universe. 

This behaviour is radically different from the corresponding local picture where 
the field, like a classical particle subject to friction, oscillates about the Eu- 
clidean/Minkowski (H — 0) minimum with progressively damped oscillations. Al- 
though the system is only mildly nonlocal (s/r rts 0.29), the increasing friction in 
the nonlocal operator eventually forbids the field to ever reach the minimum. No- 
tably, in the absence of friction the nonlocal Lorentzian solution is unbounded, as 
it undergoes wild oscillations which begin past the minimum. 33 ' 47 ' 48 ' 49 ' 500 

The branch t > corresponds to a contracting universe with scale factor a(t) = 
i" 1 / 3 , evolving towards a big crunch singularity. In this case the scalar field climbs 
the potential from its local minimum to its local maximum, pulled by antifriction. 

Since we have not solved the full cosmological inverse problem, we are not able 
to assess the physical relevance of this model as far the early history of the uni- 
verse (inflationary or pre-inflationary) is concerned. However, the first two fast-roll 

parameters are e = -H/H 2 = -3 and r? = -$/(#$) = 3[1 - (nt) 2 / '(2r)], and it 
would be hard to identify the scalar field as an inflaton generating an almost scale- 
invariant spectrum. At any rate, this may be regarded as a mathematical example 
of a bounded OSFT tachyon solution on a curved background. 



b This is not a big rip singularity 43 ' 44 ' 45 ' 46 since it happens at the infinite future where the scale 
factor is infinite. 

c The possibility that the Hubble friction be able to damp the wild oscillations was advanced in 
Refs. 4 ' 31 and realized in Refs. 14 > 50 . It would be interesting to clarify the relation between the 
solutions of Refs. 14 > 50 and ours. 
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*(U) 



V(*) 




Fig. 2. Nonlocal rolling solution (upper panel) and its potential (lower panel) for (n, a) = (2, 1). 
Here r = 1/4 and a = 1. 



3. Exponential Potential 

With little effort, one can obtain a cosmological nonlocal model with an exponential 
potential. Taking the limit a — > Eq. (1) for n = 1 + 1/ a yields 

l r 



□r(0,«^) = □e sL T(0,2 2 ) 



,-r(0„O _ 1 



(15) 



where T(a, z 2 ) — T(a) — 7(0?, z 2 ) is the upper incomplete gamma function, k 



lim Q ^o K 



c7bm/2 



1.33, and 7em is the Euler-Mascheroni constant. Equation 



(15) is satisfied with the same accuracy as the previous case; 36 therefore we refer 
again to Fig. 1. 

The potential is V(r) = — (r + e~ v )/r. Let t < 0. The universe contracts with 
scale factor a — (— t) 1 / 3 , starting at T(t — —00) = at the top of the potential 
V(Q) = -1/r. Then the field rolls down towards T(t = 0) = +00, V(+oo) = -00 
(see Fig. 3). The branch t > represents an expanding universe rolling up to the 
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Fig. 3. Nonlocal solution (top panel) for an exponential potential (bottom panel). Here r = 1/4. 

maximum of the scalar potential i.e., relaxing to the minimum in the Lorentzian 
case. 

To conclude, we have described the dynamics of two models of nonlocal scalar 
fields on cosmological backgrounds. The scalar with cubic potential drives a phase 
of superacceleration and may be identified with the bosonic OSFT tachyon. The 
cosmological friction protects the field from the wild oscillations of its Minkowski 
counterpart. 
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